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ABSTRACT 

The collisionless accretion shock at the outer boundary of a galaxy cluster should primarily heat 
the ions instead of electrons since they carry most of the kinetic energy of the infalling gas. Near 
the accretion shock, the density of the intracluster medium is very low and the Coulomb collisional 
timescale is longer than the accretion timescale. Electrons and ions may not achieve equipartition 
in these regions. Numerical simulations have shown that the Sunyaev-Zel'dovich observables (e.g., 
the integrated Comptonization parameter Y) for relaxed clusters can be biased by a few percent. 
The Y versus mass relation can be biased if non-equipartition effects are not properly taken into 
account. Using a set of hydrodynamical simulations we have developed, we have calculated three 
potential systematic biases in the Y versus mass relations introduced by non-equipartition effects 
during the cross-calibration or self-calibration when using the galaxy cluster abundance technique 
to constrain cosmological parameters. We then use a semi-analytic technique to estimate the non- 
equipartition effects on the distribution functions of Y {Y functions) determined from the extended 
Press-Schechter theory. Depending on the calibration method, we find that non-equipartition effects 
can induce systematic biases on the Y functions, and the values of the cosmological parameters ^Im, erg, 
and the dark energy equation of state parameter w can be biased by a few percent. In particular, 
non-equipartition effects can introduce an apparent evolution in w of a few percent in all of the 
systematic cases we considered. Techniques are suggested to take into account the non-equipartition 
effect empirically when using the cluster abundance technique to study precision cosmology. We 
conclude that systematic uncertainties in the Y versus mass relation of even a few percent can introduce 
a comparable level of biases in cosmological parameter measurements. 

Subject headings: cosmic microwave background — cosmological parameters — galaxies: clusters: 
general — hydrodynamics — intergalactic medium — large-scale structure of uni- 
verse 



1. INTRODUCTION 

Observational and theoretical studies have shown that 
galaxy clusters can be used as cosmological probes. In 
particular, the evolution of the galaxy cluster abundance, 
or the mass function, is sensitive to cosmological pa- 
rameters including the average matter density ^Im, the 
normalization of the power spectrum of the initial den- 
sity fluctuations and the dark energy equation of 
state parameter w (jMantz et al.l 120081 : IVikhlinin et al.l 
120091) . Except for gravitational lensing which is dif- 
ficult to do for large sample of galaxy clusters, the 
masses of galaxy clusters cannot be directly measured, 
and hence the mass function cannot be measured easily. 
Very often, the masses are estimated using mass proxies 
such as the Sunyaev-Zel'dovich (SZ) temperature distor- 
tion, the X-ray flux, or galaxy dynamics. The mass- 
observable relations have to be calibrated empirically 
or semi-empirically with numerical simulations (e.g., 
IVikhhnin et all [20091 ) . Hence, measuring cosmological 
parameters using the galaxy cluster abundance requires 
a full understanding of mass-observable relations. Even 
if the mass function and the mass-observabl e relation can 
be fit t ed simultaneously ("self-calibration" ; iLevine et al.l 
[200l lmll[200l IMaiumdar fc Mohd [200l . the correct 
form of the mass-observable relation is needed. 

A recent X-ray survey has shown that even a sam- 

|kwwong@virginia.edu, sarazin@virginia.edu, drw2x@virginia.ecru] 



pie of only 85 X-ray clusters is sufficient to provide very 
tight constrains on some cosmological parameters. For 
example, as can be measured down to 1% level in statis- 
tical uncertainty using the cluster abundance technique 
alone by assuming a flat universe with fixed dark energy 
equation of state parameter an d a prior on the Hubble 
constant fVik hlinin et al.ll2009[ ). However, the statisti- 
cal uncertainties on some other cosmological parameters 
(e.g., r^A//) are still slightly larger than 10%. Ongoing and 
future SZ surveys will detect thousands of clusters (e.g., 
Birkinshawl 119991 : iCarlstrom et al.l 120021 : iBartlett et aLl 



20081) ■ and this will significantly improve the constraints 
on cosmological parameters. Therefore, it is important 
to control the systematic uncertainties of galaxy cluster 
physics at even a percentage level. 

Because of the very long Coulomb collisional timescale 
in the low-density outer regions of galaxy clusters, it has 
been pointed out th at electron and ions there can be 
in no n-equipartition ()Fox fc Loeblll997l : lEttori fc FabianI 
Il998f ). Numerical simulations have shown that the SZ 
observables, e.g., integrated Comptonization parameter 
y), for relaxed cluster s can be biased by a few percent 
Wong fc Sarazinl[2009f ). and can potent ially be biased 
up to ^ 10% in major merging clusters (iRudd fc Nagail 
|2009| ). Specifically, the non-equipartition effect reduces 
the electron pressure compared to equipartition mod- 
els, and hence the integrated Comptonization parame- 
ter of the non-equipartition model is smaller than that 
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of the equipartition model. Although the uncertainties 
are still large, recent X-ray observations suggest that the 
electron pressure in cluster outer regions may be lower 
than that predict ed by numerical simulations assum- 
ing equipartition (iBasu et alJ |2009^ 'Geo rge et aD 120091 : 
iHoshino et 311120101 ) . Recent observations of the sec- 
ondary cosmic microwave background anisotropics with 
the South Pole Telescope (SPT) and the WMAP 7 year 
data also suggest that the electron pressure is smaller 
than the v alue predicted by hydrodynamic simulations 
(jLueker et al. 2009; Komatsu ct al. 2010). These obser- 
vational signatures are consistent with electrons and ions 
in non-equipartition, although it is also possible that the 
hydrodynamic simulations may simply overestimate the 
gas pressure. Another possibility is that heat conduction 
outside the clusters may be reducing the gas pressure 
(|Loebll2002D . 

In our previous paper, we have shown that the non- 
equipartition effect can introduce biases in the inte- 
grated SZ effect, and the biases depend on cluster 
mass and evolve with redshift in the ACDM cosmology. 
The non-equipartition model was discussed in detail in 
IWong & Sarazin (2009). In this paper, we study the im- 
pact of non-equipartition effect on precision cosmology 
studies using the non-equipartition models we have de- 
veloped. We consider only the non-equipartition effects 
associated with an accretion shock at the outer edge of 
a cluster. No significant coUisionless electron heating at 
the accretion shocks is assumed, and the intergalactic 
gas outside the coUisionless accretion shock is taken to be 
cold. These assumptions maximize the non-equipartition 
effects of the accretion shock. On the other hand, cluster 
mergers are not considered in our work. Mergers may 
increase the non-equipartition effect by a few percent, 
and hence our calculations of the biases in cosmological 
parameter estimation may still underestimate these ef- 
fects. However, the non-equiparti tion effect induced by 
mergers lasts for only 0.5-1 Gyr (jRudd fc Nagail [20091) 
which is comparable to the timescale that mergers can 
temporarily enhance the inte grated SZ effect by boost- 
ing the overall temperature (|Wik et al.l 12008'. hereafter 
WSR) . Such transient phenomenon will mainly introduce 
scatter in the Y versus mass relation, and the effect on 
cosmological parameter estimation is small in general. 
The non-equipartition effect induced by mergers may 
partially cancel out the merger boost in SZ effect, and 
hence the merger effect on cosmological parameter esti- 
mation may even be smaller. On the other hand, the non- 
equipartition effects in the accretion shock regions can 
systematically bias the Y versus mass relation for all clus- 
ters, as long as clusters are continuously accreting mate- 
rials from the surrounding which is believed to be gener- 
ally true. Moreover, systematic uncertainties in precision 
cosmology using galaxy clusters can be minimized by re- 
stricting the sample of clusters to the highest degree of 
dynamical relaxation, and hence considering the system- 
atic effects on rela xed clusters al o ne is particularly im- 
port ant. W e follow IRandall et al.l (|2002l hereafter RSR) 
and IWSRI closely to quantify the biases in cosmological 
parameter estimation using semi-analytical techniques. 
Specifically, we study the non-equipartition effect on the 
Y versus mass relation (Section (2]). Such a biased Y 
versus mass relation will affect the number of clusters 
with Y observed in SZ surveys, i.e., the Y function (Sec- 



tion In this work, the Y function is calculated using 
the e xtended Press-Schechter theory (jPress fc Schechten 
Il974f) . We consider three cases which may potentially in- 
troduce biases in the Y versus mass relations if the non- 
equipartition effect is not properly taken into account 
during the cross-calibration or self-calibration processes 
whe n us ing the galaxy cluster abundance technique (Sec- 
tion |4T]). We quantify and discuss the impact on cosmo- 
logical parameter estimation from SZ surveys by fitting 
the mass function with the biased Y versus mass relations 
(Section |4?2|) . Section [5] gives the discussion and conclu- 
sions. Throughout the paper, we assume the Hubble 
constant Hq = 71.9 hu g km s~^ Mpc~^ with /171.9 = 1. 

2. SZ VERSUS MASS CORRELATION 

The SZ effect can be characterized as the Gomptoniza- 
tion parameter, y, which is given by 



y 



P,dl. 



(1) 



where ctt is the Thomson scattering cross section, Ue is 
the electron number density, Tg is the electron tempera- 
ture, Pe = neksTe is the electron pressure, and I is the 
distance along the line of sight. The integrated Comp- 
tonization parameter, Y, is defined as the integral of the 
Comptonization parameter in equation ([1} over the area 
of the cluster on the sky 



ydA, 



(2) 



where dA is the angular diameter distance to the cluster, 
n is the solid angle of the cluster on the sky, and A is the 
projected surface area. In this paper, Y is integrated over 
the projected surface area of the cluster out to the shock 
radius. This quantity is useful for spatially unresolved 
clusters with SZ observations where the beam area covers 
the whole cluster. 

It has been shown that the integrated Comptoniza- 
tion p arameter displays a t ight correlation with cluster 
mass (jReid fc Sperge]|l2006l ). Such a tight correction is 
useful for precision cosmology, and hence a correct un- 
derstanding of the integrated Gomptonization parameter 
is important. A detailed discussion of the us e of SZ sur- 
veys to study cosmology can be found in Garls trom et al.l 
([2002). In this paper, we assume the SZ effect versus 
mass relation for the equipartition model to be the same 
as the equilibrium Y-M relation used in IWSRI . Specifi- 
cally, the equipartition SZ effect versus mass relation we 
assume is of the form 



Y^q^Nx^pix] /ifi^gMpc 



(3) 



where x = M200 / {hji t^lO^^ Mq) , N is the normaliza- 
tion constant, a is the power-law index, and p[x] is a 
13 degree polynomial in x. Equation ^ is fitted to 
the nu merica l solutions for the equilibrium Y-M rela- 
tion in IWSRI . The integrated SZ bias introduced by the 
non-equipartition effect, Y nnn-Bn/yRn versus M at differ- 
ent redshifts, is taken from 'Won g fc SarazinI (j2009( ). The 
non-equipartition SZ effect versus mass relation we used 
in this paper is hence given by 



Y„. 
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Fig. 1. — Top panel: Integrated Comptonization parameter Y 
versus mass at z = for the equipartition model used in this paper. 
Bottom panel: Integrated Y bias, (ynon-cq/yoq)wSi versus mass at 
z = and 2. 

where the subscripts "WSR" and "WS" here indicate 
that the terms are taken fro m different models in WSR 
and IWong fc SarazinI (|2009f ). respectively. In this work, 
since we are interested in the relative effects on SZ sur- 
veys and cosmological parameter estimation introduced 
by the non-equipartition effects instead of the precise Y- 
M relation which depends on details of numerical simula- 
tions, in principle, we can take any equilibrium Y-M re- 
lation from simulations and apply our non-equipartition 
bias to the equilibrium Y-M relation. The r eason for 
using the equilibrium Y-M relation in IWSRI to model 
the non-equi partition Y-M re l ation instead of the Y- 
M relation in i Wong fc SarazinI ([2009') is that the former 
relation takes into account the dependence of Y on gas 
fraction /ga., where /gas oc M^/^ for M^m > "^^^^ ^c^- 
The numerical solutions in iWong &: SarazinI (|2009D as- 
sume a constant /gas, and a self-similar argument shows 
that Y oc M^/Jj^ oc M^/^. The equilibrium Y- 
M relation in IWSRI has a power-law index close to 
2. On the other hand, the non-equipartition effects on 
the integrated Y depend weakly on /gas and hence we 
can assume the F bias of the constant /gas models in 
IWong fc Sar^ (12001) can be applied to the varying 
/gas model in |WSR|. Anoth er adv antage of using the 
equilibrium Y-M relation in IWSRI is that we can com- 
pare the effect of non-equipartition on cosmological pa- 
rameters estimations to the merger effects calculated in 
IWSRI . The equipartition Y-M relation and the inte- 
grated SZ bias introduced by the non-equipartition effect 
we used in this paper are plotted in Figure |TJ Clusters 
with higher masses are hotter, and hence, the equiparti- 
tion timescales are longer. Thus, the non- equipartition 
effects are stronger in m ore massive clusters (|Fox fc Loebl 
ll997trU^ng fc SarazinI 12^009). For our non-equipartition 
model in the ACDM universe, the integrated SZ bias 
decreases as redshift decreases. This is probably due 
to the decreasing rate of accretion onto clusters in the 
ACDM universe during the cosmological acceleration, 
whic h results in a longer tim e for electron-ion equilibra- 
tion (jWong fc Sarazinll2009D . 
In order to quantify the effect of non-equipartition on 



the SZ versus mass relation, we fit a power-law function 
to the integrated SZ bias of the form 

(rno„-cq/req)wS = AN (j^tISq j 

to all the clusters with M200 between 10^^ and 4 x 
IO^^Mq. The range is consistent with the cluster mass 
range used to study the ma, ss function in X-ray observa- 
tions (jVikhlinin et al.ll2009t ). We also fit over the wider 
range of 10^'^ to IO^^Mq for comparison. The fitted co- 
efficients AiV and Aa correspond to the biases in the 
fitted SZ versus mass relation if the non-equipartition ef- 
fect is not taken into account. The clusters in our sample 
are distributed roughly uniformly in the logarithm of the 
mass. The fitted results are listed in Table |TJ 

The effect of non-equipartition is to decrease Y for the 
high-mass clusters, and hence the effect on the SZ ver- 
sus mass relation is to lower the power-law index by 0.01 
to 0.016, which corresponds to a decrease of 0.5%~0.8% 
for the SZ versus mass relation with power-law index 
a = 2. This is comparable to the 0.8% uncertainty of 
the power-law index derived fr om simulations com bined 
with observations within M500 (|Arnaud et al.ll2009l ) . The 
normalization is biased to compensate for the change in 
power-law index, and the bias is about 3% for the scaled 
mass of 10^^ Af©. For SZ surveys which measure Y out to 
the shock radius, e.g., observations with a spatial resolu- 
tion poorer than typical shock radii, we have shown that 
the non-equipartition effect can introduce a small devi- 
ation in the measured Y-M relation. The deviations 
are small, but future SZ surveys with sufficient statis- 
tics may be able to detect such signatures. On the other 
hand, if clusters are spatially resolved and Y are mea- 
sured within i?200 , we have shown that non-equipartition 
effect is smaller than 1% for all cluste rs, and hence the 
bias in the Y-M relation is negligible (|Wong fc SarazinI 
12009; . A bias in measurement out to the shock radii but 
not within -R200 will indicate that cluster outer regions 
may be in non-equipartition. 

3. EFFECTS OF NON-EQUIPARTITION ON SZ SURVEYS 

The number of galaxy clusters expected to be found per 
unit comoving volume depends sensitively on cosmology. 
The quantity which is convenient to describe the clus- 
ter number density is the mass function, n{M, z), where 
n{M, z)dM gives the number of clusters per unit comov- 
ing volume with masses in the range M ^ M + dM, 
and z is the redshift. While the exact form of the mass 
function can be found most accuratel y from cosmolog- 
ical simulations (iSpringel et al.l 120051) , a semi-analytic 
form of the mass functi on given by the extended Press- 
Schechter theorv (Press fc Schechteil 119741: iBond et all 
IT991I: iLacev fc Cold I1993D provides a more convenient 
way to understand the dependence of the mass function 
on cosmological parameters, especially when we are in- 
terested in the relative effect instead of the precise values 
of the mass function itself. Although the Press-Schechter 
theory cannot reproduce the mass function found in cos- 
mological si mulations at very high r edshifts and low c lus- 
ter masses (jSheth fc TormenliigglfLukic et al.ll2007l ). it 
is more than sufficient over the redshifts {z = to 2) and 
masses {M = 10^"' to lO^^M©) of interest here. 
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The mass function given by the extended Press- 



Schechter theory can be written as (jPress fc Schechteil 
[1971 



npsiM,z)dM 



p Sc{z) 


da{M) 






M a^{M) 


dM 


exp 


2a2(M) 



dM . 



(6) 



where p is the current mean of the total mass density of 
the universe, a{M) is the current rms density fluctuation 
within a sphere of mean mass M, and Sc{z) is the critical 
linear overdensity required for a region to collapse at red- 
shift z. Unless otherwise specified, th e para met ers us ed 
in this paper are the same as those in lRSRl and lWSl . 

Once the Y versus mass relation is known, the distri- 
bution function of Y is given by the Y function, 

nps(r,z) =nps(M,z)^, (7) 

where nps{Y, z)dY gives the number of clusters per unit 
comoving volume at redshift z which have the integrated 
SZ parameters in the range Y Y+dY. The Y function 
for the non-equipartition models and the equipartition 
models can be related by 

dYcq 

^PS,non-cq(^^on-cq, ^) — ^PS,cq(^q, z) , (8) 

^ non-cq 

where the subscripts "eq" and "non-eq" denote the 
equipartition and the non-equipartition models, respec- 
tively. 

Figure [5] shows the Y functions for the equipartition 
and the non-equipartition models and also their ratios 
at different redshifts for the standard ACDM cosmology. 
The theoretical Y functions can be biased strongly for 
large Y and high redshift clusters. In practice, whether 
the bias can affect the observed Y functions depends on 
the number of clusters that can be observed, and this 
depends on cosmology. The maximum number of clusters 
that can be observed with Y values in the range Y — 
Y — > Y + dY and with redshifts in the range z = z — s> z + 
dz is nps(y, z) dYdV, where dV is the comoving volume 
of the universe between redshifts z and z + dz. For each 
of the four redshifts z — 0, 0.5, 1.0, and 2.0, we selected 
a redshift interval zi to z^ as given in the second and 
third columns of Table [2l Then, we defined values of Yq 
and nps,o such that J-^^ J^^ nps{Y, z) dY {dV/dz) dz = 1 
and npsfl{z) = nps{Yo,z) for the equipartition models. 
These values are listed in Table [2] and are plotted in 
Figure[2]as solid dots. For example, between z — 0.25 
0.75 for the assumed cosmology, the expected number 
of clusters with F > 6 x 10"'' /iff^g Mpc^ (M > 2 x 
10^^ h^igMQ) is one. On the other hand, the expected 
number of clusters with F = (1 ^ 2) x 10^"* ^^ng Mpc^ 
[M = (0.9 1.3) X 10^^ /ifi^gMg] is about 200 within 
the same redshifts interval, and the bias in nps is about 
5%. 



4. EFFECTS OF NON-EQUIPARTITION ON 
COSMOLOGICAL PARAMETER ESTIMATION FROM SZ 
SURVEYS 




Y (h^^ ^ Mpc') 

Fig. 2. — Top panel: Integrated Y functions of the non- 
equipartition (solid lines) and the equipartition (dotted lines) mod- 
els in the standard ACDM universe. Solid dots indicate the val- 
ues of Yq such that f^" n-psiY, z) dY {dV/ dz) dz = 1 for the 
equipartition model. Bottom panel: Ratios between the non- 
equipartition and equipartition Y functions. 

In this section, we follow a procedure similar to that 
outlined in RSl3 and IWSRI to address the impact of non- 
equipartition effects on cosmological parameter estima- 
tion from SZ surveys. Readers who a re interested in the 
technical details should refer to lRSRi and lWSRi . We out- 
line the fitting procedure used in this work and address 
the difference between the previous works below. 

As discussed in Section [3l the mass function of galaxy 
clusters is sensitive to cosmology, and hence measur- 
ing the galaxy cluster abundance at different redshifts 
can provide constraints to cosmological parameters. In 
particular, the properti es of the dark ener gy can po- 
tentially be determined (jHaiman et al.|[2"001f) . However, 
the masses of galaxy clusters cannot be directly deter- 
mined, and a mass proxy must be observed to deter- 
mine the mass through the mass-observable relation. Ex- 
amples of mass proxies are the SZ temperature distor- 
tion, the X-ray flux, and the weak leasing shear. If the 
mass-observable relations can be calibrated, this can pro- 
vide very tight constraints on cosmologic al parameters 
(Mantz_etal. 2008; Vikhli nin et all l2009| ). The mass- 
observable relations can be calibrated by numerical sim- 
ulations and/or cross-calibrations with some other ob- 
servables, and these are subjected to systematic uncer- 
tainties due to cluster physics and/or observational con- 
straints. On the other hand, the mass-observable rela- 
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tions can be simultaneously fitted with the mass function 
of galaxy clusters, and this is called "self-calibration " 
(|Levine et al.l[200a iHullMol IMaiumdar fc Mohdl200l . 
The sensitivity in constraining cosmological parameters 
and the mass-observable relations by "self-calibration" 
depends on both the form of the mass-observable rela- 
tions and the mass function of galaxy clusters. In this 
work, we are mainly interested in the biases introduced 
by the non-equipartition effect, and the effects on the 
mass-observable relations. Hence, we assume there is no 
other systematic uncertainties in the mass function in 
any given cosmology. 

4.1. Systematic Uncertainties Introduced by 
Non-Equipartition 

To address the impact of non-equipartition effects on 
cosmological parameter estimation from SZ surveys, we 
quantify the impact as biases in the cosmological pa- 
rameter estimates if the calibration of the Y versus 
mass relation does not include the non-equipartition 
effect. We generate the integrated Y function with 
the non-equipartition effects included under an as- 
sumed cosmology, and call the generated Y function 
the non-equipartition Y function. We then fit the non- 
equipartition Y function with the incorrectly calibrated 

Y versus mass relation. We consider three different cases 
for the incorrectly calibrated Y versus mass relation. For 
the first case {Case 1), we assume the Y versus mass re- 
lation is calibrated with incorrect numerical simulations 
that assume equipartition, and this incorrectly calibrated 

Y versus mass relation is used to fit the mass function. 
An example of this systematic bias might occur if the Y 
versus mass relation was extrapolated to the shock ra- 
dius using observations within a smaller radius together 
with numerical simulations assuming equipartition. In 
this case, the integrated SZ bias is simply given by 

bi = f , (9) 

where the right-hand side is the same term in equa- 
tion gl). 

For the second case ( Case 2) , we assume the Y versus 
mass relation is self-calibrated by fitting the Y versus 
mass relation and the mass function simultaneously, but 
with an incorrect functional form in the Y versus mass 
relation. We assume the incorrect functional form of the 

Y versus mass relation to be a power law in mass. In 
this case, we assume the integrated SZ bias is given by 

fe2=f%^) /f%^) , (10) 

\ ^cq / WS / \ ^cq / pifit 

where the term with the subscript "plfit" is the best-fit 
power-law relation given in equation 

For the third case {Case 3), we assume the Y versus 
mass relation is calibrated correctly at z = 0, but the 

Y versus mass relation at higher redshifts is incorrectly 
calibrated by extrapolating the calibration from that at 
z = 0. In this case, we assume the integrated SZ bias is 
given by 

63 ^ fY;^2I^\ j ^^^^ 
\ ^cq / ws,z / V ^cq / ws,2=o 



Studying the impact on cosmological parameter estima- 
tion for all of these cases is equivalently to fitting the non- 
equipartition SZ luminosity function by the equipartition 
SZ luminosity function in equation ([8|), but replacing the 
dFcq/rf^non-oq by 61,62, and &3 in equations (|9|)- (fTT|) . 

4.2. Fitting Procedures and Results 

For each case of the systematic bias we studied, we 
generate the Y function with the non-equipartition bias 
included as given in equation (|8]) by assuming the stan- 
dard ACDM cosmological model with Q,m = 0.258, 
fjg = 0.796, and a constant dark energy equation of state 
parameter w = —1. The Y function generated can be di- 
rectly calculated from the analytic Press-Schechter mass 
function in equation ([6]) and the bias ed Y ve rsus mass re- 
lations. This is simpler than those in lRSRi or lWSRi where 
merger trees were needed to generate the mass function 
in order to follow the merger history. In our work, all 
the Y versus mass relations are biased regardless of the 
merger history. We then fit the generated Y function 
with an equipartition model and find the best-fit values 
of the cosmological parameters. Clusters of galaxies can 
be used to constrain the dark energy equation of state 
parameter, w. Even the evolution of w can potentially 
be constrained. In this work, we study the constraint 
on w{z) using the form w ^ wq + wizj (1 + z)^; t he de- 
tailed explanation of this choice can be found in WSl. 
We consider three cases when fitting the cosmological pa- 
rameters: i) fitting the JIm and erg but fixing w{z) = — 1; 
ii) fitting the JIm, cs, and assuming w{z) = wq, where wq 
is a constant to be fitted; and iii) fitting the fin, erg, and 
w = Wq + wiz/{l + z)^, where wq and wi are constants 
to be fitted. The Y functions are simultaneously fitted 
at four different redshifts {z = 0,0.5,1.0, and 2.0) to 
break the degeneracy in the fitted cosmological parame- 
ters. We choose only to fit Y between 5 x 10~^ h^i.g Mpc^ 
and 5 x 10~^ /lyj^g Mpc^. The lower limit is chosen be- 
cause clusters are likely to be confused in the SZ sur- 
veys for < lOi"* /ifi^gMo HHolder et al.ll2007[) . Increas- 
ing the lower limit will make the biases in cosmological 
parameters larger since the non-equipartition effect in- 
creases with cluster mass. The upper limit corresponds 
roughly to the most massive cluster that can be formed 
in the ACDM universe. The limits are als o consistent 
with SZ surveys being done or planned fe.g.. lMelin et afl 
I2005D . We also limit the fits to values of F < Yq and 
z such that the maximum number of observable clus- 
ters J^^ nps{Y, z) dY {dV/dz) dz > 1 for most red- 
shift bins (Tabled]). 

Non-equipartition Y functions at several redshifts and 
fitted models are shown in Figure |3l The deviations in 
the non-equipartition and the fitted Y functions are small 
and only visible in the residual plots. In Case 1 and Case 
2 with the dark energy equation of state parameter frozen 
sX w — —\, for z < 1, the non-equipartition Y functions 
are higher than the best-fitted Y functions for low-mass 
clusters, and the opposite is true for high-mass clusters. 
At z = 2, the non-equipartition Y function is slightly 
higher than the best-fitted one in Case 1 with w frozen 
at w = —1, but the opposite is true in Case 2. The 
residual is similar in appearance if we free the constant 
value of w {wi frozen at zero) or allow w to vary with 
redshift in Case 1 and Case 2. In Case 3, the residual is 
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more complicated. The residuals are of the order of a few 
percent, and this will affect the estimated cosmological 
parameters as discussed below. 

The fitted results of the cosmological parameters for 
different cases we considered are summarized in Table [H 
In general, for all three systematic uncertainty cases we 
studied, when freezing w — the deviations of ^Im and 
(Tg from the assumed cosmology are < 1%. The best- 
fitted parameters happen to be consistent with the as- 
sumed cosmology, but we can see from the residual plots 
that there are clear systematic deviations by a few per- 
cent (e.g.. Figure 131). Such systematic deviations might 
be confused with the effect of non-Gaussian initial condi- 
tions on galaxy cluster mass functions, which also show 
similarly shaped s ystematic deviations (e.g.. Figure 1 in 
iFedeh et al.l[2nn9l ). 

If we free the constant value oi w = wq, non- 
equipartition effects can be significant depending on how 
the calibration is done. For Case 1, flM, erg, and wq 
now deviate by -1-3.3%, —1.6%, and —5.4%, respectively. 
This shows that ignoring the non-equipartition effects in 
the Y versus mass relation when cross-calibrating with 
numerical simulations can introduce significant biases in 
cosmological parameter estimations when one is trying to 
constrain the dark energy equation of state. Either self- 
calibrating using a power-law form in the Y-M relation 
(Case 2) or calibrating the Y-M relation correctly at 
low redshift {Case 3) can significantly reduce the biases 
in Q.M, cs, or wq (down to < 1%). 

If we allow w to vary with redshift, non-equipartition 
effects can introduce significant biases in cosmological 
parameter estimates (up to ~ 10%). Self-calibrating us- 
ing a power-law form in the Y-M relation {Case 2) can 
reduce the biases in ^Im and cg down to < 1%, but non- 
equipartition effects can still introduce a ^ 4% bias in 
the constant normalization of the dark energy equation 
of state parameter (wq) and introduce an apparent evo- 
lution of the same order. Calibrating the Y-M relation 
correctly at low redshift ( Case 3) can further reduce the 
bias in to < 1%, but again, there is still an apparent 
evolution of ^ 3%. These results show that the mass 
function technique is very sensitive to dark energy, and 
a full understanding of the systematic uncertainties in 
galaxy cluster physics is essential to constrain the dark 
energy equation of state using the mass function tech- 
nique. 

5. DISCUSSION AND CONCLUSIONS 

Numerical simulations have shown that the SZ ob- 
servables (e.g., the integrated Comptonization parameter 
Y) for rela xed clusters can be biased by a few percent 
(jWong fc Sarazin 2009), a nd potentially up to ^ 10% 
m major merging clusters dRudd fc NagaillBol ). These 
results are consistent with the SPT and the WMAP 
7 year data which indicate that electron pressures are 
smaller than predi c ted by hydrodynami c simulations 
(jLueker et al.l l2iD09l: iKomatsu et all [2010l ). A few X- 
ray observations in the cluster outer regions also show 
that the electron pressure is lower than the pressure pre- 
dictc d by numerical simulations which assume equiparti- 
tion (iBasu et al.ll2009t iGeorge et al]|2009t iHoshino et all 
120101) . However, it is also possible that numerical simu- 
lations simply overestimate the gas pressure, or that the 
gas is supported in part by other forces such as turbulent 



or cosmic ray pressure. The Y-M relation can be biased 
if non-equipartition effects are not properly taken into 
account. Precision cosmological studies using the evolu- 
tion of the galaxy cluster abundance rely on a full un- 
derstanding of the mass-observable relations, if the mass 
of the galaxy clusters cannot be directly measured as is 
usually the case in practice. We have studied systemat- 
ically the impact of biased Y-M relation introduced by 
the non-equipartition effect on SZ surveys and on preci- 
sion cosmological studies. 

While previous studies show that the Y-M relation is 
stable to complicated physical processes such as merg- 
ers a nd the power- law index in the Y-M relation is ro- 
bust (IPooleet al.ll2^007l:IWSR^ . we have shown that non- 
equipartition effect can introduce a deviation from the 
power law in the Y-M relation. We have fitted a power- 
law to the non-equipartition Y-M relation, and we found 
that there is a ^ 1% bias in the power-law index, which 
is comparable to the uncertainty of the power-law in- 
dex derived fr om simulations com bined with observations 
within M500 ()Arnaud et al.|[2009() . Such a small system- 
atic bias and deviation from the power-law form have 
important implications for SZ surveys and precision cos- 
mological studies using the SZ surveys. 

Using the analytic extended Press-Schechter theory 
to quantify the mass function of galaxy clusters, we 
have studied the non-equipartition effects on SZ sur- 
veys. We found that the Y functions can be biased 
strongly for large Y and high-redshift clusters. For ex- 
ample, the expected number of clusters with K = (1 — >■ 
2) X IQ-^h^lg Mpc2 [M = (0.9 ^ 1.3) x lO^^ /ifi^gMg] 
between z = 0.25 — J> 0.75 can be biased by ^ 5%. The 
net effect is that ignoring non-equipartition effects under- 
estimates the abundance of high-mass and high-redshift 
clusters. Cosmological parameters measured by using the 
cluster counting technique from SZ surveys will be biased 
if non-equipartition effects are not taken into account. 

We have quantified the impact of non-equipartition 
effects on cosmological parameter estimations from SZ 
surveys by the galaxy cluster abundance technique us- 
ing biased Y-M relations. We considered three poten- 
tial systematic biases in the Y-M relations if the non- 
equipartition effect is not properly taken into account 
during the cross-calibration or self-calibration when us- 
ing the galaxy cluster abundance technique. The best-fit 
cosmological parameters, f2j\f , as, and also the dark en- 
ergy equation of state parameters [w = wq+wiz/ (l + z)'^] 
using the biased Y-M relations were determined. For all 
the three methods of calibrating the F-mass relation we 
have studied, if the dark energy equation of state param- 
eter is frozen at w = — 1, we find that the best-fit and 
as are consistent with the assumed cosmology to within 
~ 1%. However, there are clear systematic deviations of 
a few percent in the fitted Y functions which may be con- 
fused wi th others effects su ch as non-Gaussian initial con- 
ditions (jFedeh et al.ll2009[ ). Models with non-Gaussian 
initial conditions predict that the actual number of clus- 
ters with higher mass can be lower than the models with 
Gaus sian initial conditions (e.g.. Figure 1 in lFedeli et al.l 
l2009f ) ; our non-equipartition model predicts an apparent 
smaller Y for high-mass clusters and this underestimates 
the number of high mass clusters if the non-equipartition 
effect is not taken into account. Note that at such small 
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Fig. 3. — Non-equipartition Y functions (crosses) and the best fitted Y functions (lines) at different redshifts are plotted on each figure 
in row 1, 3, and 5 for Case 1, 2, and 3, respectively. Residuals in the log of the non-equipartition Y functions are plotted under the 
corresponding figures. Figures in Column 1 correspond to models with dark energy parameters frozen at = 1. Figures in Column 2 
correspond to models with w fitted as constant parameters {wi frozen at zero). Figures in Column 3 correspond to models with w allowed 
to vary with redshift. Redshifts at z = 0,0.5,1, and 2 are shown in solid, dash-dotted, dashed, and dotted lines, respectively. 
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levels of systematic deviations, other systematic uncer- 
tainties such as the use of different mass functi ons (e.g., 
IPress fc SchechteJITOTa [Sheth fc Tormenlll999D may in- 
troduc e larger systematic biases. However. iFedeli et al.l 
([2009') also used the Press-Schechter mass function (as 
we do) to determine the effect of non-Gaussianity. Thus, 
we can directly compare these effects with the results of 
non-equipartition, free of biases introduced by the choice 
of mass function, and estimate the bias non-equipartition 
will introduce in efforts to detect non-Gaussian fluctua- 
tions with clusters. 

If w is fitted as a constant parameter [wi frozen at 
zero), then depending on the calibration methods, non- 
equipartition effect can introduce a few percent biases 
on the measured cosmological parameters {Case 1). Ei- 
ther self-calibrating the Y-M relation using a power-law 
form ( Case 2) or calibrating the Y-M relation correctly 
at low redshift ( Case 3) can significant reduce the biases 
in 51^/, CT8, or wq to < 1%. If we allow w to vary with 
redshift, the non-equipartition effect can introduce a bias 
in cosmological parameter of up to ~ 10% {Case 1). In 
particular, non-equipartition effects can introduce an ap- 
parent evolution in w of a few percent in all of the cases 
we considered. 

Using the cluster abundance technique alone, an X-ray 
survey with 85 X-ray clusters has already constrained 
some cosmological parameters down to 1% level in sta- 
tistical uncertainty (Vikhlinin ct al. 2009). Ongoing and 
future SZ surveys will detec t thousands of clusters (e.g., 
Birkinshawl 119991: ICarlstrom et al.l 120021: IBartlett etahl 
2008[) , and this will significantly improve the constraints 
on cosmological parameters. Therefore, it is important 
to control the systematic uncertainties of galaxy clus- 
ter physics at even a percentage level. Hydrodynamic 
simulations assuming equipartition suggest that the in- 
tegrated y is a robust mass proxy even when galaxy 
clusters are in the process of merging, and hence the 
integrated Y is ta ken to be a nearly ideal probe for cos- 
mological studies (jPoole et al.ll2007i: IWSR] ). 

Our results show that if the non-equipartition effect 
is not properly taken into account, cosmological param- 
eters can be biased significantly (up to ~ 10%). In 
order to take the non-equipartition effect into account 
when using cluster abundance to study precision cos- 
mology, the ultimate solution is to include the non- 
equipartition effect in cosmological simulations assum- 
ing the non-equipartition physics is known accurately. 
If higher resolution is needed, another approach is to 
correct the non-equipa rtition effect by perform ing ideal- 
ized simulations (e.g., iWong fc SarazinI I2009D or to re- 
simulate representative clusters taken from cosmological 
simulations including the non-equipartition effect with 
realistic assumptions. For the latter case, the non- 
equipartition effect can be taken into account together 
with other physical processes (e.g., gas depletion pro- 
cesses during the formation) which may also affect the 
Y versus mass relation. However, the above calibration 
methods by numerical simulations rely on the assump- 
tion that the non-equipartition physics is known accu- 
rately, which is in fact not the case at present. One 
of the key systematic uncertainties is the electron heat- 
ing efficiency at the collisionless shock, /3. One way to 
constrain the non-equipartition physics is to make di- 
rect observations of accretion shocks, which is currently 



not feasible. We may constrain non-equipartition physics 
based on observations of other astrophysical shocks such 
as mergers shocks and supernova remnants. However, we 
have to assume these results apply to cluster accretion 
shocks, which may or may not be the case. Another route 
might be to perform plasma simulations (e.g., particle- 
in-cell simulations) to constrain the shock physics. How- 
ever, to apply the plasma simulation results to cluster 
accretion shocks, a detailed knowledge of the pre-shock 
physics such as the magnetic field structure might be 
needed. Clearly, all of the above calculations are nec- 
essary to determine the range of the systematic uncer- 
tainties and the effects of the non-equipartition physics, 
and to constrain the form of the Y versus mass rela- 
tion. These should be studied in the near future. Until 
numerical simulations can directly determine the effects 
of non-equipartition on the Y versus mass relation from 
first principles, we suggest either to self-calibrate the Y 
versus mass relation using a power-law form at each red- 
shift bin {Case 2), or to calibrate the Y versus mass 
relation correctly at low redshift {Case 3). These will 
reduce the biases due to non-equipartition on fijv/, fs, or 
wq to better than 1%. Important biases introduced by 
other physical processes can be corrected in addition to 
the non-equipartition correction. However, if one would 
try to constrain the evolution in w to better than 1%, to- 
gether with the self-calibration method, the constraints 
from numerical simulations with uncertainties less than 
a percent level might be necessary. 

We have shown that using the cluster abundance to 
constraint the dark energy equation of state requires 
a full understanding of the systematic uncertainties in 
galaxy cluster physics. Even though we are only con- 
sidering the systematic uncertainties introduced by the 
non-equipartition effect, our results also suggest that sys- 
tematic uncertainties in the Y-M relation introduced by 
other physics of even a few percent can introduce a com- 
parable level of biases in cosmological parameter mea- 
surements. Future cluster surveys aiming to constrain 
departures from general relativity will need to control 
systematic unce rtainties down to a sub-percentage level 
(jSchmidt et al.ll2009.' ) , and hence cluster physics must be 
understood in a comparable accuracy. Future theoretical 
calculations and numerical simulations should pay par- 
ticular attention to the effects of non-thermal physics on 
the electron pressure profiles. Potential systematic un- 
certainties include conduction, turbulent pressure, mag- 
netic pressure, and relativistic pressure supported by cos- 
mic rays. Deep observations should also be carried out 
to constrain all these effects in detail for individual clus- 
ters. The outer regions of galaxy clusters are ideal sites 
for study non-thermal physics. These studies not only 
can increase our understanding of cosmology, but also 
can provide information on the physics of galaxy clus- 
ters and plasma physics under extreme conditions. 
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TABLE 1 

Effects of Non-Equipartition on the SZ-Mass 
Relation 



2 


(1-40) xlO"Af0 


(0.1-100) X lO^-* Mq 


AN 


Aa 


AN 


Aa 


0.0 


0.977 


-0.0126 


0.967 


-0.0117 


0.5 


0.971 


-0.0164 


0.963 


-0.0122 


1.0 


0.971 


-0.0152 


0.964 


-0.0109 


2.0 


0.970 


-0.0158 


0.969 


-0.0085 



TABLE 2 

Characteristic Numbers of Clusters 



z 






Yo (h-l^ Mpc2) 


Mo(hfi';glOi5MQ) 


"PS,0 Mpc 5) 


0.0 


0.00 


0.25 


6.10 X 10"'' 


2.27 


4.81 X 10-*^ 


0.5 


0.25 


0.75 


5.95 X 10-'' 


2.24 


6.01 X 10-8 


1.0 


0.75 


1.25 


1.97 X 10-'' 


1.26 


4.99 X 10-8 


2.0 


1.25 


2.75 


6.04 X lO-"' 


0.68 


2.85 X 10-^° 



Note. — The values of Yq and nps,o are defined such that 
J^^ nps{Y, z) dY (dV/dz) dz = 1 and nps,o(z) = nps{Yo,z) for the 
equipartition models for different redshift intervals zi to Zu- The value 
of Mo corresponding to Yo is related by equation 



TABLE 3 

Best-Fit Cosmological Parameters and Biases for Different Y Versus Mass 

Calibrations 



Calibration 


Qm 


0"8 


wo 


Wl 


Case 1 


0.2548(-1.2%) 
0.2665(-l-3.3%) 
0.2680(-l-3.9%) 


0.7950(-0.1%) 
0.7830(-1.6%) 
0.7811(-1.9%) 


[-1] 

-0.9464(-5.4%) 
-0.9031(-9.7%) 


[0] 
[0] 

-0.2362(+5.9%) 


Case 2 


0.2579(0%) 
0.2602(-f0.9%) 
0.2610(-fl.2%) 


0.7976(-f0.2%) 
0.7951(-0.1%) 
0.7940(-0.3%) 


-0.9890(-l.l%) 
-0.9577(-4.2%) 


[o| 

[0] 

-0.1725(-l-4.3%) 


Case 3 


0.2585(-l-0.2%) 
0.2604(-|-0.9%) 
0.2601(-l-0.2%) 


0.7953(-0.1%) 
0.7932(-0.4%) 
0.7938(-0.3%) 


[-11 

-0.9910(-0.9%) 
-1.0123(-l-1.2%) 


[o| 

[0] 

0.1167(-2.9%) 



Note. — The assumed correct cosmological parameters are JIm = 0.258, 
erg = 0.796, Wo ~ —1, and wi = 0. The bracketed values are the frozen 
values in the fits. The values in parentheses in Columns 2-4 are the percentage 
deviations of the fitted cosmological parameters from the assumed parameters. 
The values in parentheses in Column 5 are the largest percentage change in 
w between the present time {z — 0) and z — 2; this change is Aw = wi/i 
assuming w{z) — wo + wiz/{l + z)^. 



